Local simulation of singlet statistics for restricted set of measurement 
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The essence of Bells theorem is that, in general, quantum statistics cannot be reproduced by a local 
hidden variable (LHV) model. This impossibility is strongly manifested when statistics collected by 
measuring certain local observables on a singlet state, violates the Bell inequality. In this work, we 
search for local POVMs with binary outcomes for which an LHV model can be constructed for a 
singlet state. We provide various subsets of observables for which an LHV model can be provided 
for singlet statistics. 

PACS numbers: 03.65.Ud 



I. INTRODUCTION 

A violation of the Bcll-CHSH inequality [1] by statis- 
tics generated from local measurements performed on an 
entangled state shared between two spatially separated 
parties certifies such quantum state as nonlocal. The sin- 
glet state of two qubits (an EPR state) exhibits maximum 
nonlocality [2] for proper choices of local observables. Al- 
though for pure entangled states the degree of nonlocality 
is in direct proportion to the entanglement content of a 
quantum state, this is, in general, not true for mixtures 
of entangled states [3, 4]. Werner first gave the counter- 
intuitive example of mixed entangled states (popularly 
known as Werner states) [5] whose statistics when sub- 
jected to projective measurements, can be generated by 
a local hidden variable (LHV) model. A similar exam- 
ple for tripartite entangled state which can be simulated 
by a local hidden variable model was first provided in a 
work by Toth and Acin [6] . A good review of research on 
hidden variable theories can be found in [7]. 

Interestingly, Toner and Bacon [8] in the year 2003, 
gave a twist to earlier studies, by providing a model for 
singlet simulation which requires only 1 cbit of commu- 
nication supplemented with local variables. Soon after, 
Cerf et al [9] showed that 1 n/-bit (single PR-Box) is 
also sufficient for singlet simulation. Motivated by these 
works, recently, another model has been provided for sin- 
glet simulation which uses (possibly) signalling resource, 
namely S p correlations, which suggests a trade off rela- 
tion between required communication and local random- 
ness in measurement results [10, 11]. Deggore et al. [12] 
could map the problem of simulating entangled states to 
distributed sampling problems. A more through review 
of simulation of entangled state statistics from commu- 
nication complexity point of view can be found in [13]. 



Few other recent works [14-17] show that lack of free 
will can also be considered as a resource for singlet sim- 
ulation. There are also some efforts in solving the dif- 
ficult problem of simulating multipartite entanglement 
and non-maximally bipartite entangled states either by 
use of communication or by nonlocal (no-signaling) re- 
sources [18, 19]. All these various approaches have been 
deepening our understanding about quantum correlation 
and its use as a physical resource in various information 
processing tasks. 

As of providing local variable models for class of en- 
tangled states, in a seminal work in the year 2002, Bar- 
rett [20] generalized the work of Werner [5], by con- 
struting a LHV model for any positive-operator-valued 
measurements at the expense of the weight associated 
with singlet in Werner state. Motivated by these works 
we, in this paper, pose the problem from opposite di- 
rection i.e. rather than weakening the (singlet) state we 
search for the class of (weakened) dichotomic observable 
(POVM) for which local model can be provided. In par- 
ticular, here we provide the subset of the most general 
two outcome measurements represented by positive op- 
erator value measure (POVM) and present local models 
for singlet statistics generated from them. We provide 
some sets of local observable which are optimal for the 
protocol we have suggested. First we show that, if ob- 
servable on any one side is sufficiently restricted (deviates 
from ideal projective measurement), resulting statistics 
for the singlet state has a local hidden variable model. 
Next, we provide another model which is symmetric in a 
sense that observable on both the sides are put to a simi- 
lar type of restriction. Finally, we identify a more general 
set of observable for which LHV models exists with some 
further restrictions. Before we derive our results, in the 
followings section, we give a mathematical description of 
a general two-outcome POVMs. 
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II. GENERAL TWO-OUTCOME POVM 

Generalized quantum observables are described by 
POVMs [21]. For finite, say n, outcome measurements 
on a d-dimensional state space a POVM is a collection 
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of selfadjoint operators {Ei} acting on a complex Hilbert 
space C d satisfying the conditions: (i) < Ei < I for all 
i, and (ii) A = ^> where i £ {1, 2, n}. A measure- 
ment of such an observable {Ei} on a quantum state p 
results in any one of the n possible outcomes; the proba- 
bility of an occurrence of z-th outcome (termed as click- 
ing of i-th effect) is Tr\pE^\. A subclass of these type 
of general measurement has an interesting physical inter- 
pretation as unsharp spin properties, introduced by P. 
Busch [22, 24]. 

In this work, we consider general two-outcome POVMs 
{E, I — E} acting on C 2 (state space of a qubit). Effect 
E is characterized by some parameters, say, ao G K 
(a scalar) and a G M 3 (a vector). We denote norm of 
a by p. Then, the selfadjoint property along with the 
condition < E < I implies that E can be expressed as 

E=~[a I + pa-a] (1) 

< a < 2 (2) 
< p < min{ao, 2 — ao} (3) 

where a-a = a x a x +a y a y +a z a z . Then, the corresponding 
operator I — E is also selfadjoint and satisfies the require- 
ment < I - E < I, Thus, the Eq.(l) along with the 
conditions (2) and (3), supplemented with an arbitrary 
direction a, completely determine a two-outcome POVM 
{E, I — E} acting on C 2 . The region feasible for param- 
eters ao and pi for defining such an effect E{ao,p, a} is 
illustrated in Fig.(l). 
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[22, 24] first showed that a subclass of general two- 
outcome POVMs can be interpreted as measurement of 
unsharp-spin property of spin-l particles. Under the con- 
dition of rotation covariance, parameters {ao, p} is decou- 
pled from a which can then be interpreted as orientation 
of the measuring device. Further, condition of symmetry 
under a rotation ir of the measuring device gives ao = 1. 
Thus, effect operators for an unsharp spin observable is 
of the form E±(a) = \[I ± pa ■ a]. The spectral decom- 
position of positive operators E±(a) is 

where g [I + a.a] and h [I — a.a] are one dimensional spin 
projection operators on the Hilbert space C 2 . Now, the 
quantity (-^) can be suitably interpreted as degree 
of reality (unsharpness) of outcomes obtained from a spin 
measurement along direction a. From this representation 
it is clear that the POVM {E$(a), Et{a)} is a smeared 
version of the projective measurement { | [/ + a.a] , | [/ — 
a.a]} — in case of projective measurements the unsharp 
parameter p — 1. 

Another important property is that under suitable con- 
ditions two POVMs can be jointly measurable [25]. Two 
POVMs of the form {E U I - E x } and {E 2 ,I - E 2 } are 
jointly measurable if there exits a four-outcome POVM 
{E12, E12, E X 2, E12} such that it can reproduce the cor- 
rect marginals, i.e., E\ = Ei 2 + E 12 - and E 2 — En +Ei 2 - 
For unsharp spin observable it has been shown that [22] 
(also see the review [23]) two observables parameterized 
by, say, (pi,di) and (p2,a 2 ) are jointly measurable if 
and only if ||/xidi + ^2^211 + Umi^i — ^2 a 2\\ < 2. On 
considering unsharp parameter for both the spin observ- 
ables to be same i.e., p\ = p 2 , along with the fact 
|| Si + 0.2 1| + ||oi — 0.2 1| < 2~\/2 for any pair of unit vec- 
tors a\ and a 2 , it is easy to conclude that if the unsharp 
parameters p\ = p 2 < then joint measurement of un- 
sharp spin property can be realized for any such pair of 
directions. 



III. LHV MODEL FOR SINGLET STATISTICS 
FOR TWO OUTCOME POVMS 



FIG. 1: Parameter ao (p) varies along the horizontal (vertical) 
axis. Any point (ao, p) laying in the shaded triangular region 
POI (together with an arbitrary parameter a) determines a 
two-outcome POVM {E,I — E}. Points on the dashed line PU 
represent unsharp spin measurements. Point P(l, 1) represent 
ideal projective measurements. 

An interesting application of general two-outcome 
measurements considered here is in the study of spin 
properties of spin-^ systems. In this context, P. Busch 



Suppose, two spatially separated parties Alice and Bob 
share one qubit each from a singlet state 

pAB = ^[1 ® I - a x ® a x - a y ® a y ~ a z ® a z ]. 

Let Alice's (Bob's) observable be a most general two- 
outcome POVM Ea[cio, pa, a] (Eslbo, pb, o]), defined by 
Eq.(l). If the effect Ea(b) clicks we denote the outcome 
by 'yes' otherwise 'no'. Then joint outcome probabilities 
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are following: 



'yes' with a probability 



P AB (yes,yes) = 


4 [ao&o 


- ha^bo- ■ b] 


P AB (yes,no) = 




- bo) + MAM-Ba ' b] 


P AB {no,yes) = 




ao)bo + ^a^bo. ■ b] 


P AB (no, no) = 


>- 


Oo)(2 - 6 ) - VAVBa ■ b] 



(4) 



A. Models for two-outcome measurements 

Violation of the Bell-CHSH inequality [1] implies that 
there can be no LHV model for the singlet statistics 
generated by projective measurements by both the par- 
ties. Therefore, the statistics of the singlet can have a 
LHV model only if general two outcome POVMs con- 
sidered here are restricted (deviate from ideal projec- 
tive measurements) in some way or the other. Fol- 
lowing Werner's local model for some mixed entan- 
gled states [5], we provide two LHV models for sin- 
glet state under certain restrictions on parameters of 
two outcome POVMs. In both type of models vectors 
A = (sin 6> cos </>, sin sin 0, cos 6*) uniformly distributed 
over the unit sphere, are the local variables preshared 
between Alice and Bob. 



1. A fully biased model M/^: 

Let, Bob's observable EB(bo, UB,b) satisfy restriction 
Ms < ^min{6o •> 2 — b } but there is no restriction on Al- 
ice's observables, see Fig(2). 

Alice, for her observable Ea [do , ha , o] > declares l yes' 
with a probability 



yAt. «o . 1 



Pi (yes) 



\la cos a 



(5) 



where a is angle between direction a and A. On the 
other hand, for observable £'b[6o, ub, b], Bob declares 



^f (yes) = - hb sgn(cos/3) 



(6) 



where (3 is the angle between the direction b and A, and 
sgn(x) = +1 (-1) for x > (x < 0). 

The joint probability of the outcome (yes, yes), can 
be calculated from 

P& B (yes, yes) = [ p(X) if (yes) if (yes) dX (7) 



where p(X) is the considered (uniform) distribution of the 
hidden variable A. Evaluating the above integral gives 

p ihv (yes, yes) = ■- [a b - haVbo. ■ b] (8) 

which exactly matches with the quantum mechanical pre- 
diction for the outcome (yes, yes). The desired quantum 
mechanical probabilities for the other possible outcomes 
easily follows, for example, P t AB (yes, no) is obtained sim- 
ply by replacement if (yes) — > P?(no) = 1 — P?(yes) 

A A A 

in the integrand of the Eq.(7). 

2. A fully symmetric model M/ s : 

Let Alice's and Bob's observables satisfy following re- 
striction (see Fig. (3)) 



il a < -^=min{a , 2 - a } 

Ms < ^=min{6 ,2- 6 } 

Alice declares 'yes', for her observable Ea[clq, HA,a], 
with a probability 



1 



(9) 



where a is angle between direction a and A. 





FIG. 2: Alice's (Bob's) parameters can take values from the 
dark gray triangular region on left (right). Alice's parameters 
ao and fj,A can take any possible value, but Bob's parameters 
bo and iib are restricted to come from the region MOI. 



FIG. 3: Alice's (Bob's) parameters can take values from the 
dark gray triangular region on left (right). Alice's parameters 
ao and pa as well as Bob's parameters bo and pb are restricted 
in the same way and come from the triangular region MOI on 
the left and right respectively. 
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Where else Bob declares 'yes' 
Eb\Pq, Hb, b], with a probability 



for his observable 



P*(yes) = ^--^MBSgn(cos/3) (10) 



where j3 is the angle between the direction b and A, and 
sgn(:r) = +1 (-1) for x > (x < 0). Like in the fully 
biased model Mfb, we find that this model (M/ s ) also 
simulates the correct statistics for the singlet. 

If we consider unsharp spin properties on both sides 
with uniform value of fiA and \j,b and also assume any 
pair are jointly measurable [22, 23] on both sides, then 
the conditions of the model Mf s are automatically sat- 
isfied and hence this LHV model M/ s can simulate the 
singlet statistics for any arbitrary pair of respective di- 
rections a and b for Alice and Bob. 



B. Measure of restriction on observable 

By considering that observables of Alice and Bob are 
picked from a uniform distribution of all possible two- 
outcome POVMs, we can define a measure r for restric- 
tion on the observables of any of the two parties in the 
following way. (see Fig. (2) and Fig. (3)) 



1 - 



Area 



Area 







x 100 



(11) 



Now, one can easily calculate that in the model Mfb 
(M/g), there is 0% (29.3%) restriction on Alice's observ- 
ables wherelse Bob's observables are restricted by 50% 
(29.3%). 

Another interesting observation is that model M f b and 
M.f s belong to a general class of LHV models, {M K : n > 
0}. 

Under the restrictions, 



fiA < t min{ao, 2 — ao} 

Ms < ^- min{6 ,2 - b } 
Zk 

Alice declares "yes" with a probability 

A Z Zk 

Bob declares 'yes' with a probability 



Pi (y^) = ~ 



bo 
2 



k/ib sgn(cos/3) 



(12) 
(13) 

(14) 

(15) 



For any nonnegativc value of n we get a LHV model — 
M fb {Mfs) correspond to K = 1 (k = -j=). In Fig(4) 
the two curves show the % restriction on Alice's and 
Bob's observable for LHV models corresponding to dif- 
ferent values of n. The intersection point of two curves 
correspond to the symmetric model Mf s . Observe that 
k = j correspond to another fully biased model, say 



M'/b, which is same as Mfb except that conditions on 
Alice's and Bob's observables are interchanged. In fact 
all the models for which n e (0, |] U [1, oo) are fully bi- 
ased models, however one can immediately observe that 
within this subclass, either Mfb or M' fb is sufficient to 
simulate any other fully biased model. Thus only the 
subclass {M K : k € [1/2, 1]} contains tight LHV models 
in a sense that they can capture any varying degree of 
restrictions on Alice's and Bob's observables. 
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FIG. 4: Black (Gray) curve show percentage restriction on 
Alice's (Bob's) observable for LHV models corresponding to 
different values of k. The intersection point of two curves at 
k = correspond to the completly symmetric LHV model 
M/ s . Fully biased model Mf b (M' f b ) correspond to k = | 
(k = 1). All the models for n G (0, |] U [1, oo) are fully 
biased. 



C. A different class of model 

In previous cases, we put restriction on the observable 
separately on both sides. Now we put the following re- 
striction on the observable where one of the restrictions 
involves parameters of both sides; 



1 „ 1 

- < a < 2 - - 
■q T) 

HaUb < 7T min{6o, 2 - b } 
2/7 



(16) 
(17) 



where 77 > 1. Alice and Bob now can simulate singlet 
statistics according to the following protocol 



Pi (yes) = 



1 

2t? 



cos a 



>b, b 



Pi (yes) 



tjha^b sgn(cos/3) 



(18) 



(19) 



But this model is obviously non-local as Bob's output 
involves parameters of observable on both sides. But 
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the model can be made local for a given rj and fixed 
[i A- In this case there is no restriction on direction 
a of Alice's POVM. It might seem that by increasing 
the value of r\ the range of clq can be extended but 
then due to the condition (17), the range of /is is also 
restricted accordingly. So in some sense in this model clq 
and hb maintain a complementary relation for a given 6 n - 



cases where LHV simulation is possible for singlet state. 
We find the optimal set of two outcomes observable for 
which singlet simulation by LHV is possible under the 
suggested protocol. It is also interesting that for uniform 
unsharp parameter, the joint measurability of unsharp 
spin property on both sides implies LHV model for sin- 
glet. It will be interesting to study whether the set can 
be enlarged with respect to different LHV model. 



IV. CONCLUSION 
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